A Trick to Improve the Efficiency of Generating Unweighted Be Events from 

BCVEGPY 
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In the present paper, we provide an addendum to improve tlie efficiency of generating unweiglited 
events witliin PYTHIA environment for tlie generator BCVEGPY2.1 [C.H. Cliang, J.X. Wang 
and X.G. Wu, Comput.Phys.Commun.174, 241(2006)]. This trick is helpful for experimental 
simulation. Moreover, the BCVEGPY output has also been improved, i.e. one Les Houches Event 
common block has been added so as to generate a standard Les Houches Event file that contains 
the information of the generated Be meson and the accompanying partons, which can be more 
conveniently used for further simulation. 
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NEW VERSION PROGRAM SUMMARY 

Title of program : BCVEGPY2.1a 

Program obtained from : CPC Program Library 

Reference to original program : BCVEGPY2.1 

Reference in CPC : Comput. Phys. Commun. 175, 
624(2006) 

Does the new version supersede the old program: No 

Computer : Any LINUX based on PC with FORTRAN 
77 or FORTRAN 90 and GNU C compiler as well 

Operating systems : LINUX 

Programming language used : FORTRAN 77/90 

Memory required to execute with typical data : About 
2.0 MB 

No. of bytes in distributed program, (including 
PYTHIA 6. 4) : About 1.5 MB. 

Distribution format : .tar.gz 

Nature of physical problem : Hadronic Production of Be 
meson and its excited states. 

Method of solution : To generate weighted and un- 
weighted Be events within PYTHIA environment 
effectively. 



Restrictions on the complexity of the problem : Hadronic 
production of (c6)-quarkonium via the gluon-gluon 
fusion mechanism are given by the 'complete calculation 
approach'. The simulation of Be events is done within 
PYTHIA environment. 

Reasons for new version : More and more data are 
accumulated at the large hadronic cohider, it would be 
possible to make precise studies on Be meson properties, 
such as its lifetime, mass spectrum and etc.. The 
BCVEGPY has been adopted by several experimental 
groups due to its high efficiency in comparison to that 
of PYTHIA. However, to generate unweighted events 
with PYTHIA inner mechanism as programmed by the 
previous verison is still time-consuming. So it would 
be helpful to improve the efficiency for generating 
unweighted events within PYTHIA. Moreover, it would 
be better to use an uniform and standard output format 
for further detector simulation. 

Typical running time : Typical running time is machine 
and user-parameters dependent. I) To generate 10^ 
weighted S'-wave (c5)-quarkonium events (IDWTUP=3), 
it will take about 40 minutes on a 1.8 GHz Intel P4- 
processor machine. II) To generate unweighted S'-wave 
(c6)-quarkonium events with PYTHIA inner structure 
(IDWTUP=1), it wiU take about 20 hour on a 1.8 GHz 
Intel P4-processor machine to generate 1000 events. HI) 
To generate 10^ unweighted S'-wave (c6)-quarkonium 
events with the present trick (IDWTUP=1), it will take 
17 hour on a 3.16 GHz Intel E8500 processor machine. 
Moreover, it can be found that the running time for the 
P-wave (c6)-quarkonium production is about two times 
longer than the case of S-wave production under the 
same conditions. 



* email: 
t email: 



xianyou.wangiaicern 



wuxgQcqu.edu.cnT 



1^ 



Keywords : Event generator; Hadronic production; Be 
meson; Unweighted events. 

Summary of the changes : 



2 



(1) The generator BCVEGPY [J-II has been pro- 
grammed to generate B^. events under PYTHIA envi- 
ronment [3], which has been frequently adopted for the- 
oretical and experimental studies, e.g. Refs. [B-flsj. It 
is found that each experimental group shall have its own 
simulation software architecture, and the users will spend 
a lot of time to write an interface so as to implement 
BCVEGPY into their own software. So it would be bet- 
ter to supply a standard output. 

The LHE format becomes a standard format p^ . 



which is proposed to store process and event information 
from the matrix-element-based generators. The users 
can pass these parton-level information to the general 
event generators like PYTHIA and HERWIG "20"] for 
further simulation. For such purpose, we add two 
common blocks in genevent.F. One common block is 
called as bcvegpy_pyupin and the other one is writejhe. 
bcvegpy_pyupin, which is similar to PYUPIN subroutine 
in PYTHIA, stores the initialization information in the 
HEPRUP common block. 



INTEGER MAXPUP 
PARAMETER (MAXPUP=100) 

INTEGER IDBMUP,PDFGUP,PDFSUP,IDWTUP,NPRUP,LPRUP 
DOUBLE PRECISION EBMUP,XSECUP,XERRUP,XMAXUP 
COMMON/HEPRUP/IDBMUP(2),EBMUP(2),PDFGUP(2),PDFSUP(2), 
&IDWTUP,NPRUP,XSECUP(MAXPUP),XERRUP(MAXPUP), 
&XMAXUP (MAXPUP) ,LPRUP (MAXPUP) 



The writeJhe, which is similar to PYUPEV subroutine in the HEPEUP common block, 
in pythia, stores the information of each separate event 



INTEGER MAXNUP 
PARAMETER (MAXNUP=500) 

INTEGER NUP,IDPRUP,IDUP,ISTUP,MOTHUP,ICOLUP 

DOUBLE PRECISION XWGTUP,SCALUP,AQEDUP,AQCDUP,PUP,VTIMUP, 

&SPINUP 

COMMON/HEPEUP/NUP,IDPRUP,XWGTUP,SCALUP,AQEDUP,AQCDUP, 
&IDUP(MAXNUP) ,ISTUP (MAXNUP) ,MOTHUP (2,MAXNUP) , 
&ICOLUP(2,MAXNUP),PUP(5,MAXNUP),VTIMUP(MAXNUP), 
&SPINUP(MAXNUP) 



The above information can be separately exported 
to two files 'bcvegpy.init' and 'bcvegpy.evnt'. Then 
with the help of the PYTHIA subroutine PYLHEF, 
we can combine these two files into a single LHE file 
'bcvegpy.lhe'. 

(2) As for previous version of BCVEGPY [IHi, its main 
concern is to improve theoretical estimations, i.e. to ob- 
tain a precise result as fast as possible. For the purpose, 
the VEGAS [2l| running together with IDWTUP=3 is 
usually appreciated. Firstly, by running VEGAS, one 
can obtain a sampling importance function, and then 
the phase-space points are generated according to the 
relative importance of their corresponding differential 
cross-sections. This sampling importance function is 
useful, since the integrand for the phase space integra- 
tion can be regularized for obtaining a high precision 
total/differential cross section [IJ]. Secondly, by set- 
ting IDWTUP=3, all parton-level events when input to 



PYTHIA shall always be accepted, i.e. almost all the cal- 
culated points are effective ones and shall be evolved as 
final events. As a combination of these two choices, one 
can generate events within the so-called importance sam- 
pling approach, i.e. the parton-level events have been ef- 
fectively input into PYTHIA with proper weighted phase 
space. This is one of the reason why BCVEGPY runs so 
fast 1. 

Experimentally, it is the unweighted inputs other than 
the weighted inputs that are usually needed. By default 
BCVEGPY generates unweighted events following the 
same scheme of PYTHIA. For the purpose, one have to 
set PYTHIA parameter IDWTUP=1, then aU events in a 



Nevertheless, the dominate reason for its quick running is due to 
the use of improved helicity amplitude approach to deal with the 
hard scattering amplitude as described in detail in Ref. [ij. 
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PT(Bc)(GeV/c) ii(Bc) P^(Bc)(GeV/c) ii(Bc) 



FIG. 1: Comparison of the normalized transverse momen- 
tum (Pt) and pseudo-rapidity (77) distributions between 
the weighted ones (IDWTUP=3) and the un-weighted ones 
(IDWTUP=1), which are represented by solid line and dot- 
ted lines respectively. 



run are treated on an equal footing. And the parton-level 
events are uniformly generated within the whole phase 
space, but shall with a weight (XWGTUP) when input to 
PYTHIA, and then these events are accepted or rejected 
with probability XWGTUP/XMAXUP (XMAXUP is 
the maximum weight of the process). As most of the 
calculated parton-level events are rejected by this crite- 
ria, such a naive scheme is usually time-consuming. For 
example, by setting IDWTUP=1, it takes about 20 hours 
on a 1.8 GHz Intel P4-processor machine to generate 1000 
Be events. It is too slow if one wants to simulate large 
amount of events. 

So, we present a new scheme to improve the efficiency 
of generating unweighted events within PYTHIA envi- 
ronment. For the purpose, we first unweight the events 
by using the standard hit-and-miss technique p2j . Then, 
we set IDWTUP=1, but we replace XMAXUP to be 
a reference weight that can be generated by earlier ex- 
ploratory run ^. Each parton-level event shall result in 
a ratio between the weight XWGTUP to the reference 
weight. When such ratio is smaller than one, the event 
shall be retained or rejected based on whether the ratio is 
greater or less than the return value of random number 
generating function PYR(O). When the ratio is larger 
than one, its integral part stands for the basic number 
the same events to be generated. And the basic number 
shall be increased by one or unchanged based on whether 
its decimal part is greater or less than a random number 
PYR(O). It is found that a larger reference weight means 
more events to be rejected and less efhciency. So, to set 
a proper reference weight is a tricky problem, which re- 
quires exploratory running under one's own experimental 
cuts and should has negligible derivation from theoreti- 
cal estimations. For example, one may observe that for 



FIG. 2: Comparison of the normalized transverse momen- 
tum (Pt) and pseudo-rapidity (rj) distributions derived by 
BCVEGPY and PYTHIA, which are represented by red-solid 
line and blue-dotted lines respectively. 



ATLAS transverse momentum cut Pxcut = 4 GeV and 
the rapidity cut ycut =2.5, only few events' weight are 
larger than the value of 10^ pb. So we can set this valve 
as the reference weight, which is reasonable for the LHC 
Be production simulation. 

The code for the newly suggested unweighted scheme 
is put in the main program 'bcvegpy.F'. As a cross 
check, we present the weighted and unweighted Be trans- 
verse momentum Px-distribution and pseudo-rapidity rj- 
distribution in Fig. ([I}, which are represented by solid 
line and dotted lines respectively. The unweighted dis- 
tributions are obtained by setting IDWTUP=1 and us- 
ing the new strategy described above, while the weighted 
distributions are obtained by setting IDWTUP=3. It 
can be found that the weighted and unweighted distri- 
butions are coincide with each other. Since the new un- 
weighted scheme is time-saving in comparison with that 
of the previous one, so it can be adopted for generating 
large amount of events, which may suit the needs of the 
newly accumulating data at LHC. 

Moreover, to compare and validate the result of 
BCVEGPY, we also use PYTHIA's inner mechanism to 
generate Be meson at the very beginning. For such pur- 
pose, we set MSEL=5 to generate bb first, and then add 
one filter (PDGID=:541) that searches and stores every 
event that shall result in at least one Be meson. Next, 
we use the external decay generator EVTGEN [2^, 
which can be running under the PYTHIA environment, 
to simulate the Be decay channel Be J/ip + n. The 
remaining parameters for PYTHIA are taken to be their 
default values. Moreover, in doing the calculation, we 
set the collision energy of LHC to be the present value 
of 7 TeV. The comparison is presented in Fig. ([2]), which 
shows that the normalized distributions for PYTHIA 
and BCVEGPY agree with each other well. Such an 
observation also agrees with the previous conclusion 
drawn by Ref . [11] . 



^ The function of the reference weight is similar to that of 
XMAXUP, and for convenience, we still use the label XMAXUP 
to store its value. 
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